Collisional decay of a strongly driven Bose-Einstein condensate 
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We study the collisional decay of a strongly driven Bose-Einstein condensate oscillating between 
two momentum modes. The resulting products of the decay are found to strongly deviate from the 
usual s-wave halo. Using a stochastically seeded classical field method we simulate the collisional 
manifold. These results are also explained by a model of colliding Bloch states. 
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The decay of many-body states coupled to a quasi- 
continuum of collisional products is a topic of great ex- 
perimental and theoretical interest Q, 0] • Experimental 
Bose-Einstein condensation (BEC) allows us to investi- 
gate this decay in detail by use of highly controlled opti- 
cal lattice potentials. Both the coherent evolution of the 
condensate in the lattice, and the nature of the quasi- 
continuum can be manipulated and quantified 

The finite lifetime of perturbative bulk excitations in 
BEC, namely Beliaev and Landau damping, has been ex- 
tensively studied and is rather well understood. 
These studies were extended recently to the ground state 
of a BEC in an optical lattice, and weak excitations over 
such a state, using band theory formulation |& Ifl. Ilflj. 

Coherent Rabi oscillations of the condensate between 
two (or more) macroscopically-populated momentum 
states can be driven by a strong moving optical lattice 
potential. These oscillations are described as beating be- 
tween Bloch states belonging to different bands of the 
lattice 0, . Due to interaction nonlinearity, a super- 
position of two such beating Bloch states no longer solves 
the time-dependent Schrodinger equation, and thus the 
expected spectra and dynamics are richer than for a sin- 
gle Bloch state 0. The decay of these excited states 
cannot be described by mean-field theory nor as an inter- 
action between perturbative Bogoliubov quasi-particles 
[l^ as in the Beliaev formalism. 

In this letter we study the collisional decay of such 
a strongly driven BEC undergoing coherent Rabi oscil- 
lations between momentum states by a resonant two- 
photon Bragg process [13 . We measure a clear deviation 
of the collisional products from the s-wave halo observed 
for collisions of a weak excitation with the BEC 1]. Us- 
ing a stochastically seeded classical field Gross-Pitaevkii 
equation (GPE) simulation we observe similar decay 
dynamics. These results are then explained by a model 
which includes collisions between Bloch states of the op- 
tical lattice as a perturbation. 



* Current address: Time and Frequency Division NIST 325 Broad- 
way Boulder, Colorado 80305, USA 



As in H3, our nearly pure (~ 95%) BEC of TV = 
1.6(±0.5) x 10 5 87 Rb atoms in the \F,m f ) = |2,2) 
ground state, is formed in a magnetic trap with radial 
and axial trapping frequencies of uj r = 2tt x 226 Hz and 
ou z = 2tt x 26.5 Hz, respectively, leading to a healing 
length £ = 0.23 jam. The condensate is driven by a pair of 
strong Bragg beams counter-propagating along the axial 
direction z, with wavenumbers kdi = — and k^i = ^l, 
with kL = 27r/780 nm. The laser frequency is red- 
detuned 44 GHz from the 87 Rb D2 transition in order 
to avoid spontaneous emission. The depth of the result- 
ing optical lattice potential is characterized by the two- 
photon Rabi frequency Qd- For strong excitations the 
mean-field shift is largely suppressed 12], due to tem- 
poral averaging of the shift to zero during a cycle of 
oscillation, and hence the frequency difference between 
the dressing beams (in the laboratory frame) is set to 
Sd = 2tt x 15 kHz, the free-particle resonance. This leads 
to Rabi-like oscillations between the momentum states 
k = and k = 2k^. 

The oscillation in the momentum of the atoms is ap- 
parent in Fig. d where we plot the measured average 
momentum per particle in the z direction as a function 
of time, extracted from time of flight images. The oscil- 
lation frequency, as obtained by a decaying sinusoidal fit, 
is tt d /27r = 8.6 kHz. 

In the strongly driven condensate, both finite size 
broadening, and inhomogeneous density broadening are 
greatly suppressed |l2| . Therefore the decay of the oscil- 
lations is mostly due to the collisions between atoms in 
momentum modes and 2hkL- The products of such a 
collision have an average momentum of hkL, and do not, 
in general, oscillate any more in momentum space. For 
a Bogoliubov excitation, which is a weak excitation of 
momentum 2hkL over a large condensate of zero momen- 
tum, the collisional products are known to be located on 
a shell in momentum space, known as the s-wave halo. 
This shell is the surface in momentum space, conserv- 
ing both energy and momentum for the collision. Due 
to the Bogoliubov dispersion relation, which is quadratic 
for 2/cl£ > 1, this shell is nearly spherical. 

In our experiment, however, the condensate is strongly 
driven at large momentum 2k = 3.9, and collisions 



2 



CM 



CD 

E 
o 

E 

CD 
D) 
03 
s— 

<D 
> 
03 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



\ 



■ * 







200 400 
time (jisec) 



600 



FIG. 1: Average momentum per particle contained in the 
atomic cloud as a function of time (in units of 2hkL, in the 
laboratory reference frame). Oscillations are due to a strong 
moving optical lattice which is suddenly switched on, leading 
to Rabi oscillations between momentum wavepackets. The 
decay in the oscillations (fit by the dashed line) is mainly due 
to collisions which deplete the condensate. The arrow marks 
the point at which Fig. E|a) (and the subsequent theoretical 
figures) are taken. 



occur within the lattice potential. Consequently, the Bo- 
goliubov description is no longer adequate. This is clearly 
visible in Fig. Efa), which shows an absorption image ob- 
tained after a resonant dressing pulse lasting t = 660 /is. 
Upon comparison with the s-wave sphere obtained when 
two condensates collide [Fig. EJb)], one sees a clear shift 
of the collided atoms towards the center of the sphere. 
To quantify this difference we employ computerized to- 
mography to extract the radial dependence of the density 
from the column-density available in the absorption im- 
age 0. In Fig. |2fc) we plot the radial distribution of 
atoms over a small slice in z. This inward shift of col- 
lided atoms is robust and is clearly observed for different 
values of fid and t. 

We qualitatively simulate this collisional decay using 
the stochastically seeded classical field method, in 2D, 
as it was developed recently for colliding BECs 
In this method the initial seed of fluctuating random 
amplitudes of the bosonic field is added to the ground 
state of the condensate in the harmonic trap. Then the 
moving lattice potential is switched on suddenly (as in 
the experiment). Matter wave mixing between the con- 
densate momentum wavepackets and the seeded quasi- 
continuum drives the collisional decay of the oscillations, 
without any need for further gross numerical interven- 
tion. Fig. OJa) shows the resulting momentum distribu- 
tions for collisions occurring within the lattice, and Fig. 
Efb) shows the collisional manifold when only a weak lat- 
tice is present. We note that the overall agreement be- 
tween the simulation and experiment takes into account 




FIG. 2: (a,b) Absorption images after t to f = 38 ms time of 
flight following a resonant dressing pulse, (a) Strong pulse 
fld/^Tr — 8.6 kHz, pulse width t = 660 fis (b) Weak pulse 
^d/27r < 2 kHz, pulse width t = 370 fis. Dotted circles 
represent the predicted s-wave shell. x r = HkL/M x ttof is 
the ballistic expansion distance of an atom with wavenumber 
kL (in lattice frame of reference). The collisional manifold 
for the strong pulse is clearly shifted inwards as compared to 
that of the weak pulse, which agrees with the expected s-wave 
shell, (c) The density distribution along the y-axis obtained 
by computerized tomography 18] of the data in (a,b) averaged 
over a slice marked by the vertical dashed lines. The solid line 
is for the strongly driven BEC of (a), and the dotted line is 
for the weakly excited BEC (b). The collisional products of 
the strongly driven BEC are clearly driven towards the center, 
while those of the weakly driven BEC are concentrated on the 
s-wave sphere (y = x r ). 



many possible systematic effects such as the harmonic 
confinement in the radial dimension and finite time of 
the Bragg pulse. We also see that the essential physics is 
qualitatively captured here, even though the simulation is 
not in 3D. The correlations between counter-propagating 
momentum wavepackets, clearly visible in the simulation, 
are not visible in the experiment due to the fact that the 
experimental absorption images integrate over an addi- 
tional dimension, making this signal difficult to observe. 
The mean-field broadening of the experimental time-of- 
flight images also leads to some additional smearing. 

To obtain an intuitive model that still captures the 
essence of these phenomena, we neglect inhomogeneous 
and finite-size effects and choose our frame of reference 
as moving with the optical lattice of the Bragg lattice 
beams at velocity v = —hk^jM along z. Our system 
is then described by the many-body time-independent 
Hamiltonian 
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FIG. 3: (a,b) Momentum distributions generated by the clas- 
sical field GPE simulation (a) Strong pulse Qd = 8.6 kHz, 
pulse width t = 660 /us (b) Weak pulse Qd < 2 kHz, pulse 
width t = 370 /xs. The strong lattice leads to a clear shift 
of the collisional products inward as compared to the weak 
pulse, and deviates strongly from the s-wave collisional sphere. 
Note that the simulation also generates momentum wavepack- 
ets with clear number correlations, (c) The density distribu- 
tion along the y-axis, of the simulation, averaged over a slice 
marked by the vertical dashed lines. The solid line is for the 
strongly driven BEC of (a), and the dotted line is for the 
weakly excited BEC (b). The inward shift, observed in the 
experiment [Fig. Efc)], of the decay products is even more 
pronounced here. 



of Bloch states \n) b = ^2 p a n ((2p + l)fcrW2p + l)k L ) w , 
which diagonalize the Hamiltonian (|2j) [ll|, where n is 
the Bloch band index. The lattice momentum hq = hkL 
remains unchanged and is therefore omitted. The sub- 
script w and b indicate whether the quantum numbers 
in the ket describe the wavenumber of a plane- wave or a 
Bloch band index. 

In the weak lattice limit, we arrive at the two state 
result \kL) w = 1/a/2 (|l) b + |2) 6 ), of a two level system 
undergoing Rabi oscillations with frequency f^. This 
two mode picture is still useful even for stronger lattices, 
since the energy separation between the lower two bands 
to the third band, on the edge of the Brillouin zone, is 
typically larger than the discussed here. Therefore, 
the higher bands are only weakly occupied by the system. 

We now consider the mixing of Bloch states due to 
atomic interactions. In order to describe collisions, we in- 
clude the interaction term in the Hamiltonian (QJ , which 
scatters two atoms from the populated states. We fo- 
cus on the processes in which atoms are scattered into 
the quasi-continuum of unpopulated states, neglecting 
scattering into populated states (forward scattering) [2l| . 
Since collisions are binary the Hilbert space is reduced to 
a two particle space. Hamiltonian JJJ can be rewritten 
in the basis of the Bloch Hamiltonian (0) as 

H= \vi;v2){vi\V2\[E Vl +E V2 \ + 

^l|k 1 >(z, 2 |k 2 )(k 3 |z,3)(k4|^> X 

ki ...k4,Z^i ...U4 

kl,^2)(^3,^4^(ki+k2)-(k 3 +k 4 )- (3) 
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where a^ak) is the creation (annihilation) operator of a 
particle with wave- vector k, and g is a constant describ- 
ing the s-wave interactions. The relative momentum in 
the experiment is sufficiently low to avoid higher par- 
tial wave collisional terms [3 E3- Neglecting for the 
moment the interaction term, Eq. (QJ simplifies into a 
one-dimensional single particle Hamiltonian, 



H 
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hVtd cos(2kLz). 



(2) 



According to Bloch's theorem, a state \k) w is only cou- 
pled to states \k + 2pkL)w, where p is an integer. In the 
moving frame of reference, the stationary BEC has mo- 
mentum HkL and is situated on the Brillouin zone bound- 
ary. The initial kinetic energy of the condensate is there- 
fore in the lattice energy gap. To calculate the conse- 
quent dynamics we span the state \kL) w by the new basis 



Here V{ stands for all quantum numbers of a Bloch state 
ni,Qi,ki±. E Vi — E n ^ q + (hki±) 2 /2M is the energy of the 
noninter acting Bloch state, where E n ^ q is an eigenvalue 
of the Bloch Hamiltonian @ and is the part of k^ 
which is perpendicular to z |22|. The inclusion of the 
quantum numbers k^ is necessary since collided atoms 
gain momentum which is not along z. We treat the col- 
lision term in Eq. © by use of perturbation theory [2^ . 
That is, we assume the system is undergoing coherent os- 
cillations in time due to the lattice potential, and study 
the perturbative collisional products that are created by 
the interaction Hamiltonian. 

The existence of two macroscopically occupied, distinct 
energy states, implies several decay routes for the colli- 
sional term, and the subsequent energy and momentum 
conservation manifold are split. Specifically, splitting 
arises from the energy difference between the case where 
both colliding atoms are from the n = 1 band and when 
both are from the n = 2 band. Due to symmetry, de- 
structive interference suppresses collisions in which there 
is initially one atom in each of the two bands. 

The prediction of our model for the momentum dis- 
tribution of the collisional products is plotted in Fig. 0] 
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FIG. 4: Column momentum distribution for the same pa- 
rameters as for Fig. EJa) and|5Jb) respectively calculated by 
our colliding Bloch state model (Eq. only the scattered 
atoms are shown). Dotted line represents the s-wave shell. 
The momentum distribution is only roughly equivalent to the 
spatial distributions of Fig. |2]due to interactions during ex- 
pansion, although the overall shape of the collisional manifold 
is reproduced. 

The presence of the optical lattice is found to drive 
the collided atoms towards the center, as expected from 
the inner shell of the splitting. The amplitude of the 
outer shell decreases rapidly as Qd increases, and is ex- 
perimentally unobservable for our parameters. In the 
experimental data additional effects such as inhomoge- 
neous broadening and mean- field expansion effects, that 
are neglected in this model, tend to broaden the shell of 
atoms, and blot out the splitting for smaller f^. 

One intriguing prediction of the model, is that the de- 
cay rate as a function of time will deviate significantly 
from that predicted by the Fermi golden rule. This can 
be explained in the time domain by the oscillatory behav- 
ior of the coherent evolution, leading to a complementary 
oscillation in the rate of production of collided atoms. 

Another important result of the model is that as we 
increase Q,^ the overall decay of the coherent evolution is 
accelerated. This is due to the large number of additional 
decay pathways that are switched on by the presence of 
a deeper lattice. 

In conclusion, we measure the collisional decay of a 
driven condensate and show that it deviates from the 
usual s-wave sphere. This result is modelled by using the 
stochastically seeded classical field method applied to the 
Gross-Pitaevskii simulation of the experiment. The main 
features of the collisional decay manifold are captured by 
a simple model, which treats the interactions as binary 
collisions between single particle Bloch states. 
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